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Executive  Summary 


This  report  originated  in  the  H°°  Research  Initiative  of  the  Office  of  Naval  Research 
and  the  In-House  Laboratory  Independent  Research  (ILIR)  Program  of  Space  and 
Naval  Warfare  Systems  Center  San  Diego  (SSC  San  Diego).  These  programs  focused 
on  H°°  engineering  for  fleet  applications — wideband  impedance  matching  and  wide¬ 
band  amplifier  optimization.  Research  in  these  applications  produced  several  papers 
[33],  [32],  [2],  [3],  four  patents,  a  book  [1],  and  sparked  the  Defense  Advanced  Research 
Projects  Agency’s  interest  in  digital  H°°  engineering. 

H°°  engineering  computes  the  best  possible  performance  bounds.  For  example,  a 
wideband  antenna  should  be  matched  to  the  line  impedance  to  minimize  power  and 
prevent  amplifier  burnout.  The  circuit  designer  matches  the  antenna  by  searching 
for  a  lossless  2-port  that  minimizes  the  Voltage  Standing  Wave  Ratio  (VSWR).  Tra¬ 
ditionally,  the  circuit  designer  guesses  a  2-port  topology  and  then  optimizes  over  its 
circuit  elements.  This  process  is  repeated  over  various  topologies,  hoping  that  a  2-port 
that  is  “good  enough”  turns  up.  In  contrast,  H°°  engineering  computes  the  smallest 
VSWR  attainable  by  any  lossless  matching  2-port  independent  of  circuit  topology 
[21],  This  best  VSWR  provides  an  absolute  benchmark  to  assess  candidate  circuits, 
and  brings  some  order  to  matching-circuit  selection.  Moreover,  the  H°°  methods  also 
produce  some  information  on  an  optimal  2-port  circuit.  If 


S 


•Si  1  <Si2 
S21  S2-2 


is  the  unknown  scattering  matrix  of  an  optimal  2-port,  then  the  H°°  methods  compute 
sn-  Consequently,  techniques  that  dilate  the  passive  1-port  sn  into  a  lossless  2-port 
solve  the  synthesis  problem  [5]. 

Dilations  are  basic  to  electrical  engineering,  signal  processing,  and  operator  the¬ 
ory.  This  report  makes  explicit  the  algebraic  structure  of  these  matrix  dilations  by 
the  Cosine-Sine  Decomposition  (CSD).  The  CSD  provides  a  unifying  computational 
framework  for  parameterizing  all  unitary  dilations  of  a  given  matrix,  parameterizing 
all  contractive  dilations,  and  generalizes  to  produce  all  J-unitary  dilations.  These 
dilations  are  foundational  for  the  synthesis  problem.  Moreover,  the  CSD  applies  to 
several  signal-processing  problems  [29],  [35],  [16].  Consequently,  the  CSD  is  a  simple 
and  flexible  technique  that  can  be  applied  to  theoretical  and  computational  dilation 
problems. 
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1  Matrix  Dilations 


Matrix  dilations  can  be  approached  in  several  ways.  The  historical  approach  could 
start  with  Halmos’  original  work  on  unitary  and  normal  dilations  [17].  The  harmonic 
analysis  approach  could  start  from  the  work  of  Nagy  &  Foia§  [13]  that  links  analytic 
functions  and  dilations.  Another  approach  follows  from  the  link  between  operators 
and  electrical  circuits  discovered  by  Helton  [19];  in  this  approach,  circuit  synthesis 
is  equivalent  to  either  unitary  [28]  or  J-unitary  matrix  dilations  [20],  [21].  These 
approaches  are  only  a  few  ways  to  access  the  massive  literature  on  dilation  theory. 
This  report  approaches  dilation  theory  using  the  Cosine-Sine  decomposition  (CSD). 
The  focus  is  on  matrix  dilations  so  the  algebraic  patterns  clearly  stand  out. 

Let  A  be  a  given  M  x  M  matrix.  The  dilation  problem  is  to  parameterize  all 
unitary  dilations  U a  of  matrix  A:  find  all  matrices  B  G  CMx7V,  C  G  <CNxM ,  and 
D  G  CNxN  so  that  the  dilation 


UA  = 


A  B 
C  D 


is  unitary: 

Im  0 
0  In 

Here  Im  denotes  the  M  x  M  identity  matrix  and  the  superscript  H  denotes  the 
Hermitian  or  conjugate  transpose.  The  CSD  parameterizes  all  unitary  dilations  of  A 
and  reveals  how  the  dilation  of  smallest  size  is  encoded  in  A. 

A  unitary  dilation  is  not  the  only  type  of  dilation  that  can  be  considered.  A 
significant  generalization  replaces  the  unitary  equality 


KUa 


Im+n 


with  the  contractive  inequality 


U a  U a  <  Im+n- 


In  the  latter  case,  Ua  is  called  a  contractive  dilation  of  A.  The  CSD  parameterizes 
all  these  contractive  dilations.  The  simple  patterns  give  a  short  proof  of  Parrot’s 
Theorem,  which  has  connections  to  H°°  theory  and  electrical  engineering.  It  is  in 
electrical  engineering  that  the  CSD  generalizes  to  its  J-unitary  version.  In  this  setting, 
the  problem  is  to  find  all  dilations  of  A  that  are  J-unitary  or  hyperbolic: 


VaJUa  =  J-= 


Im  0 
0  —In 
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Electrical  engineers  routinely  use  a  natural  mapping  between  unitary  matrices  and 
hyperbolic  matrices.  Under  this  mapping,  the  CSD  naturally  turns  into  the  Hyper¬ 
bolic  Cosh-Sinh  decomposition  (HCSD).  From  this  HCSD,  all  J-unitary  dilations  are 
obtained.  Thus,  the  CSD  admits  a  hyperbolic  generalization  with  engineering  appli¬ 
cations.  We  conclude  by  pointing  out  that  several  other  matrix  decompositions  also 
generalize  to  hyperbolic  versions  with  signal-processing  applications. 


2  Notation 


The  set  of  all  complex- valued  M  x  N  matrices  is  denoted  by  CMxiV.  if  x  G  CMxN . 
X  is  said  to  have  size  M  x  N ,  which  will  be  denoted  X  ~  M  x  N.  The  Hermitian  or 
conjugate  transpose  of  X  is  denoted  by  XH .  The  group  of  M  x  M  unitary  matrices 
is  denoted  by 

U{M)  :=  { U  G  CMxM  :  UHU  =  IM}- 
A  diagonal  matrix  0  is  denoted  by 


0  =  diag(0i,  •  •  • ,  9m ) 


‘  0 
0  02 


0 

0 


0  0  ...  9m 


With  a  slight  abuse  of  notation,  cos(@)  denotes  the  diagonal  matrix 


cos(0i)  0 
,  0  cos(02)  . 

cos(0)  = 

I 

0  0 

Similarly,  sin(@)  denotes  the  diagonal  matrix 

’  sin(0i)  0 

°  sin(02) 

sm(0)  = 

0  0 


0 

0 

cos(0 


M) 


0 

0 


sin(0 


M 


3  Cosine-Sine  Decomposition 


If  a  matrix  A  is  square  (N  =  M),  one  unitary  dilation  of  A  is  the  Halmos  dilation 
[17,  Chapter  23]: 


A  ~(IM  -  AAHf/2  ' 

(. Im-AhAY /2  AH  ’ 
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where  the  positive  semidefinite  square  root  is  selected.  Halmos  points  out  that  this 
dilation  has  a  nice  geometric  interpretation.  If  A  has  a  singular-value  decomposition 
A  =  U  cos(0)VH,  where  U  and  V  are  unitary  matrices,  then  Ua  generalizes  the 
classic  rotation  matrix  as  follows: 

_  U  0  cos(0)  —  sin(0)  V  0  H 

A  0  V  sin(0)  cos(0)  0  U 

This  is  a  special  case  of  the  CSD. 


Theorem  1  (CSD)  [34,  page  37],  [15,  page  77]  Let  the  unitary  matrix  W  G  U(M  + 
N )  be  partitioned  as 


W 


w12 

c  M 

w21 

w22 

CN 

CM 

CN 

If  N  >  M,  then  there  are  unitary  matrices  Uu,  Vn  G  U{M )  and  unitary  matrices 
U 22,  V22  G  U(N)  such  that 


where  C  >  0  and  S  >  0  are  diagonal  matrices  satisfying  C2  +  S2  =  Im- 


Thus,  the  CSD  simultaneously  produces  the  singular-value  decompositions  (SVDs) 
for  Wn,  W21,  ITi 2 ,  and  IT22  from  the  sine  and  cosine  matrices.  The  converse  param¬ 
eterizes  the  unitary  dilations  of  A. 


4  Parameterizing  Unitary  Dilations 

Given  the  SVD  of  a  matrix  A,  the  CSD  shows  us  how  to  get  a  unitary  dilation  of  A. 
The  utility  of  the  CSD  is  that  all  unitary  dilations  of  A  with  N  >  M  can  be  obtained 
this  way. 

Corollary  1  [23,  Problem  1.6.21]  Let  A  G  CMxM  be  a  contraction.  Select  any 
singular-value  decomposition 

A  =  Uu  cos(@)Dq) 


3 


where  0  =  diag(#i,  •  •  • ,  6m)  is  sorted  in  ascending  order:  0  <  6\  <  ■  ■  ■  <  9m  <  vr/2. 
If  N  >  M,  then  all  unitary  dilations  of  A  with  D  e  £NxN  are  parameterized  as 


A  B 
C  D 


Uu  0 
0  U 22 


cos(0) 

—  sin(@)  0 

sin(0) 

cos(0)  0 

0 

O 

'  Yn 

0  ' 

0 

y22  _ 

where  U2 1,  Y22  G  U(N). 


Proof:  If  Ua  has  the  given  CSD,  then  Ua  is  a  unitary  dilation  of  A  and  demonstrates 
that  unitary  dilations  of  A  exist  for  N  >  M.  Conversely,  if  Ua  is  a  unitary  dilation 
of  A  with  N  >  M,  the  CSD  gives 


UA  = 


Xu  0 

0  X22 


cos(@) 

—  sin(@)  0 

sin(0) 

cos(0)  0 

0 

l 

0 

'  Yn 

0  ' 

0 

y22  _ 

lH 


where  Jn,  Yn  G  U(M),  X22,  Y22  G  U(N),  and  A  =  Xn  cos(Q)Yff .  The  unicity 
of  the  SVD  allows  us  to  set  Xn  =  Uu  and  Yu  =  Vn  and  use  only  U22  and  V22  to 
parameterize  the  unitary  dilations.  To  demonstrate  this  claim,  observe  that  A  has 
the  SVDs 

A  =  Un  cos (0)Cn  =  Xn  cos(0)X1f . 


So,  how  unique  are  Uu  and  Vu  compared  to  Xn  and  Yn?  Write 


cos(0) 


c  1 1  !f 


C2-/77 


CrIt, 


where  C\  >  c2  >  ■  ■  ■  >  Ck  >  0  are  the  distinct  non-zero  singular  values  of  multiplicity 
rrik-  By  [23,  Theorem  3.1.1'],  there  are  unitary  matrices  Wi,  W2,  ■ . . ,  IT Y ,  Ei,  and 
E2  such  that 


'  TYi 

'  W1 

W2 

w2 

F  = 

and  G  = 

WK 

WK 

Ei  . 

e2  _ 

link  the  two  SVDs  of  A  as  Xn  =  UnF  and  Yn  =  Yn G.  By  construction,  both  F 
and  G  commute  with  cos(@)  and  sin(0),  and  satisfy 

cos(@)  =  Fcos(Q)Gh  =  Gh  cos (0)F. 


4 


Substitution  gives 


UA 


where 


UuF  0 
0  X22 


t/n  0 
0  X22 


t/n  0 
0  X22 


Cn  0 
0  U22 


cos(0) 

—  sin(O) 

0 

sin(O) 

0 

cos(O) 

0 

§ 

O  1 

GHVR  0 
0  Y» 


'  F  cos(Q)GH 

— F  sin(0) 

0 

sin(0)G" 

cos(0) 

0 

L  0 

0 

1 

cos(0) 

L 

—  sin(0)F 

0 

GH  sin(0) 

GH  cos (0)F 

0 

0 

0 

1 

[  YR 

0  ' 

0 

YH 
I22  J 

r  vH 
V11 

0  ' 

0 

yh 

I22  J 

cos(0) 

—  sin(0) 

0 

sin(0) 

0 

cos(0) 

0 

0 

In-m 

vR  o 

o  VR 


U'22  ~  X22 


GH  0  ' 
0  In-m 


V22  =  Y, 


22 


Fh 

0 


0 

In-m 


Because  X22,  Y22  0  U(N)  are  arbitrary,  U22  and  V22  are  also  arbitrary.  Thus,  we  may 
fix  Uu  and  Vn  to  parameterize  Ua  using  only  U22  and  V22  0  U(N).  /// 


The  proof  shows  how  the  non-uniqueness  in  the  SVD  of  A  (the  matrices  F  and 
G)  can  be  peeled  off  and  then  cast  into  the  arbitrary  U2 2  and  V22  matrices.  Thus,  U\\ 
and  Vn  may  be  fixed  so  that  the  dilations  are  parameterized  only  by  U22,  V22  0  U(N). 

For  dilations  with  N  <  M,  consider  the  following  numerical  example.  Partition 
the  random  unitary  matrix  Ua  as  follows: 


'  -0.3944 

-0.2202 

0.1686 

0.0281 

-0.5633 

0.4974 

0.4495  ' 

r 

-1 

-0.5276 

-0.0880 

-0.5915 

0.1517 

0.1399 

0.3175 

-0.4697 

A 

B 

0.6000 

-0.0909 

-0.2244 

0.6823 

-0.2040 

0.2708 

-0.0318 

C 

D 

— 

-0.0484 

0.1497 

-0.3991 

-0.0607 

-0.7027 

-0.5589 

-0.0777 

0.3307 

-0.7570 

-0.2090 

-0.5057 

-0.0567 

0.0589 

-0.1069 

L 

-0.2331 

0.2001 

-0.4216 

0.4028 

-0.2224 

-0.5647 

0.3707 

-0.3375 

0.3144 

0.1602 

-0.4437 

0.2563 

0.5342 

0.5230 

Corollary  1  would  dilate  A  to  a  10  x  10  unitary  matrix.  However,  A  has  the  singular- 
value  matrix 


h  0 

0  cos(@) 


1.0000 

1.0000 

1.0000 

0.8771 

0.2452 

The  two  smallest  singular  values  clue  us  that  the  5x5  matrix  A  came  from  the  7x7 
unitary  matrix  Ua-  More  generally,  the  number  of  singular  values  that  are  strictly 
less  than  1  actually  encode  the  size  of  the  smallest  unitary  dilation. 


Corollary  2  [23,  Problem  1.6.21]  Let  A  E  CMxM  be  a  contraction  with  SVD 


A 


Uu 


II  0 

0  cos(@) 


VH 


11  > 


5 


where  0  =  diag(6b,  •  •  • ,  Ok)  and  0  <  9i  <  ■  ■  ■  9k  <  7r/2.  If  N  <  M,  then  all  unitary 
dilations  of  A  with  D  ~  N  x  N  must  have  N  >  K  and  take  the  form 


'  A 

B  ‘ 

'  Uu 

0  ‘ 

C 

D 

0 

U22  _ 

Im-n 

0 

0 

0 

cos(T) 

—  sin(T) 

0 

sin(T) 

cos(T) 

Vn  0 

0  V22 


where  U2 2,  V22  €  U{N)  and  T  = 


°N-K  0 

0  © 


Proof:  The  CSD  permits  us  to  select  the  decomposition 


D  C 
B  A 


X22  0 

0  Xn 


cos(T) 

sin(T) 

0 

—  sin(T) 

cos(T) 

0 

0 

0 

Im-n 

'  e22 

0  ' 

0 

H 


where  T  =  diag^i  0  <  r/q  <  . . .  <  <  7t/2;  X22,  Y22  G  U{N)\  Xn, 

Yu  e  U{M).  By  some  non-trivial  relabeling,  we  may  write 


"  A 

B  ' 

'  Xn 

0  ' 

C 

D 

0 

X22 

Im-n 

0 

0 

0 

cos(\f') 

—  sin(’I') 

0 

sin(’I') 

cos(4') 

Yu  0 
0  Y22 


(1) 


By  assumption,  A  has  K  singular  values  strictly  less  than  1  and  L  singular  values 
exactly  equal  to  1.  Then  K  +  L  —  M  and  M  —  N  <  L.  Then  all  dilations  must  have 
N  >  M  —  L  =  K.  The  unicity  of  the  singular  values  [23,  page  146]  and  the  ordering 
of  the  'fn's  permit  us  to  write  T  as  stated.  When  N  =  K,  it  follows  that  T  =  0 
so  that  the  smallest  dilations  exist.  Equation  1  now  permits  the  application  of  the 
unicity  arguments  from  Corollary  1.  Thus,  we  may  fix  Xn  and  Yu  as  the  unitary 
matrices  Uu  and  Vu  from  the  SVD  of  A  and  sweep  out  all  dilations  of  A  by  sweeping 
over  X22,  Y22  e  U(N).  /// 


The  important  special  case  of  Corollary  2  is  the  parameterization  of  the  unitary 
dilations  of  smallest  size: 


'  A 

B  ‘ 

'  Un 

0  ‘ 

C 

D 

0 

U22  _ 

II 

0 

0 

0 

cos(@) 

—  sin(0) 

0 

sin(@) 

cos(0) 

C11 

0 


0 


H 


V22 


where  |cos(6lm)|  <  1.  This  case  explains  the  preceding  numerical  example  where 
A  ~  5  x  5  had  Jl  ~  3x3  and  0  ~  2  x  2.  Thus,  A  has  unitary  dilations  of  size  7x7 
or  larger. 
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5  Parameterizing  Contractive  Dilations 


Given  a  matrix  A  that  is  a  contraction,  the  problem  is  to  hnd  all  dilations 


Ta  = 


A  B 
C  D 


that  are  also  contractions:  TffTA  <  Im+n ■  The  CSD  provides  a  straight-forward 
parameterization  of  all  the  TVs  by  compressing  a  unitary  dilation. 


Corollary  3  (Adapted  from  [7],  [30],  [27],  [36],  [If,  Corollary  3.5])  Let  A  e  CMxM 
be  a  contraction.  Select  any  singular-value  decomposition 

A  =  U\\  cos(Q)Vf[, 


where  0  =  diag(#i,  •  •  • ,  9m)  is  sorted  in  ascending  order:  0  <  9i  <  ■  ■  ■  <  9m  <  tt/2. 
Then  all  contractive  dilations  of  A  with  D  ~  N  x  N  may  be  parameterized  as 

cos(0)  —  sin(@)  0 

sin(@)  cos(0)  0 

0  0  I2N 


A 

C 


B 

D 


Un 

0 


0 

U22 


0 

U23 


r  vH 

vn 

0 

0 

VM 

0 

\ 

where  [U22  U23],  [ V22  V23 ]  G  P+M ( M  +  2 N).  Here,  Pn  denotes  the  orthogonal 
projection  onto  the  first  N  components  of  CM+2N . 


Proof:  Let  Pm+n  denote  the  orthogonal  projection  onto  the  first  M  +  N  components 
of  £2M+2N '  By  Corollary  1,  A  has  unitary  dilations  of  the  form 


'  Uu 

0 

0 

cos(0) 

—  sin(0) 

0 

r  vH 

vn 

0 

UA  = 

0 

U22 

U23 

sin(0) 

cos(0) 

0 

0 

VH 

v22 

VH 

17 .32 

0 

U32 

U33 

0 

0 

I2N 

0 

VH 

v23 

VH 

v33  J 

Then  TA  =  Pm+nU a\Cm+n  is  a  contractive  dilation  of  A  with  D  G  CNxN .  Conversely, 
let  T a  be  a  contractive  dilation.  By  Corollary  1,  the  unitary  Halmos  dilation 


Uta 


Ta  -(Im+n  -  TAT f)1/2 

(Im+n  ~  T»Ta)1/2  Tf 


admits  the  factorization 


'  Uu 

0 

0 

cos(0) 

—  sin(0) 

0 

r  vH 

vn 

0 

Uta  = 

0 

U22 

U23 

sin(0) 

cos(0) 

0 

0 

vs 

VH 

v32 

0 

U32 

U33 

0 

0 

I2N 

0 

V2H3 

VH 
v33  J 

Then  Ta  =  Pm+nUa\Cm+n .  /// 


This  whole  apparatus  generalizes  to  operators  on  Hilbert  spaces  [14].  A  good 
representative  of  the  CSD  in  action  on  a  Hilbert  space  is  Parrot’s  Theorem. 
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Theorem  2  (Parrot)  [7],  [30],  [27],  [f]  LetH  mid  K.  be  Hilbert  spaces  with  orthogo¬ 
nal  decompositions  H  =  7i\  ©  7d2  and  1C  —  JC\  ©  /C2.  Let  A  :  7i\  — >  LC\,  B  :  H2  —■ ►  1C\, 
C  :  7i\  — >  JC2,  be  fixed  operators.  For  D  :  7d2  — >  1C 2,  define  To  :  7i  — >  1C  as  the 
operator 

1C 1 

/c2  . 

Hi  h2 


A  B 
C  D 


Then 


inf{||TD||}  =  max 


A  0 
C  0 


A 

0 


(2) 


The  CSD  provides  a  proof  of  Parrot’s  Theorem  for  the  matrix  case.  Without 
loss  of  generality,  scale  the  right  side  of  Equation  2  to  1  so  that  ||TD||  >  1  for  any 
operator  D.  Equality  is  demonstrated  by  finding  a  contractive  To-  The  scaling  also 
forces  both  matrices  in  the  right  side  of  Equation  2  to  be  contractions.  By  Corollary  3, 
both  matrices  admit  the  representations: 


'  A 

0  ' 

'  £7ncos(0)^f 

0  ' 

C 

0 

_  t722  sin  ( 0 )  Vfii 

0 

'  A 

B  ' 

'  Uu  008(0)^  — C/nsin(0)H2f  1 

0 

0 

- 1 

O 

O 

_ 1 

Combining  both  representations  gives  the  contraction 


A  B 
C  D 


Unc  os(0)V(f  -t/nsin(0)C2f 
U22  sin(0)Vlf  f/22cos(0)C2f 


and  proves  Parrot’s  Theorem. 

Meinguet  [27]  offers  a  fine  exposition  of  Parrot’s  Theorem  and  its  applications.  A 
key  application  is  Nehari’s  Theorem  [36],  [30]  that  is  the  foundation  of  H°°  engineering 
[22],  A  link  with  dilation  theory  can  be  traced  as  follows.  By  1970,  Nagy  &  Foia§ 
[13]  developed  a  dilation  theory  for  analytic  functions  whose  values  are  contractions 
on  a  Hilbert  space.  In  1972,  J.  W.  Helton  [19]  connected  their  dilation  theory  to 
the  main  realizability  theorem  of  electrical  engineering.  By  1982,  Helton  [21]  and  his 
colleagues  had  made  deep  connections  between  operator  theory,  electrical  engineering, 
and  control  theory.  These  applications  enriched  the  operator  theory  with  significant 
generalizations.  In  particular,  the  dilation  theory  admits  a  nice  generalization  to 
dilating  J-unitary  or  hyperbolic  matrices.  Electrical  engineering  problems  provide  an 
exellcnt  motivatation  to  consider  the  hyperbolic  matrices  and  J-unitary  dilations. 


6  Hyperbolic  Matrices  in  Electrical  Engineering 


A  basic  object  in  electrical  engineering  is  the  iV-port.  The  iV-port  is  a  “black  box” 
with  N  pairs  of  wires  sticking  out  of  it.  The  word  “port”  means  that  each  pair  of 
wires  obeys  a  conservation  of  current — the  current  flowing  into  one  wire  of  the  pair 
equals  the  current  flowing  out  of  the  other  wire.  Figure  1  shows  a  2-port  with  voltage 
sources  driving  each  port.  The  AT-port  is  the  collection  of  voltage  v  and  current  i 
vectors  that  can  appear  on  its  ports  [20] .  Thus,  the  iV-port  is  really  a  subset  Af  of  a 
larger  voltage-current  space,  typically  an  L 2  space. 


2 -Port 

s= 

Ul  H2 
s2  \  S22 

o- 

+ 


o 


Figure  1:  2-port  voltages,  currents,  and  scattering  parameters. 


Physical  assumptions  about  the  Y-port  translate  into  geometric  statements  about 
A f.  For  example,  a  linear  iV-port  is  equivalent  to  the  subset  Af  being  a  linear  subspace. 
Under  more  restrictive  assumptions,  the  AT-port  can  be  the  graph  of  a  linear  operator. 
For  example,  if  the  Y-port  that  relates  voltage  and  current  as  v  =  Z\  is  characterized 
by  the  impedance  matrix  Z\ 


Likewise,  if  i  =  Yv,  the  Y-port  is  characterized  by  its  admittance  matrix  Y : 


Af  = 


V 

Yv 


Not  all  Y-ports  necessarily  admit  impedance  or  admittance  matrices — classic  exam¬ 
ples  are  open  circuits,  short  circuits,  and  transformers.  However,  if  each  port  is 
connected  to  a  voltage  source  and  series  resistor  rn,  the  claim  is  that  any  linear,  time- 
invariant,  solvable  Y-port  is  characterized  by  its  N  x  N  scattering  matrix  S  [5],  [6], 
[20].  Specialized  to  the  2-port  in  Figure  1,  the  scattering  matrix  maps 


1*1 

Sll 

Sl2 

CL\ 

<N 

rO 

_ 1 

S21 

S22 

d2 
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where  the  incident  signal 


a  =  (i?g  1/2v  +  i?g/2i)/2 


and  the  reflected  signal 

b  =  (-Ro  1/2v  —  i?o/2i)/2 

are  computed  from  the  voltage  v  and  current  i  via  the  normalizing  matrix 


Rq  — 


r'i  0 

0  r2 


No  loss  of  generality  is  incurred  by  taking  Ro  as  the  identity  matrix:  R$  =  R. 

The  power  P  consumed  by  the  2-port  is  provided  by  Balabanian  and  Bickart  [6, 
pages  241-242]: 

P  =  9ft  [vHi]  =  ||a||2  —  ||b||2  =  a  H  (I2  —  SH  S)  a.  (3) 


If  the  2-port  consumes  no  power  (P  =  0)  for  all  its  voltage  and  current  pairs,  the 
2-port  is  lossless.  By  Equation  3,  a  2-port  is  lossless  if  and  only  if  its  scattering 
matrix  S  is  unitary:  SH S  =  J2.  How  do  dilations  fit  into  circuit  theory?  If  Port  2 
is  terminated  in  resistor  r2,  then  the  reflectance  Si  looking  into  Port  1  is  Si  =  Sn. 
The  circuit  synthesis  problem  is  the  converse:  given  the  reflectance  si,  End  all  lossless 
2-ports 


S 


Sl 1  Si2 

S  2 1  s22 


with  si  =  Sn . 


Thus,  circuit  synthesis  is  a  problem  in  dilation  theory  [21],  [20],  [28]. 


Closely  related  to  the  scattering  matrix  is  the  chain  scattering  matrix  0  [18,  page 
148]: 


'  bl  ' 

\  On 

6 12 

a2 

=  0 

02  1 

CL\ 

(N 

_ i 

O22 

.  &2 

1 

<N 

rO 

_ 1 

When  multiple  2-ports  are  connected  in  a  chain  as  in  Figure  2,  the  chain  scattering 
matrix  of  the  chain  is  product  of  the  individual  chain  scattering  matrices. 


When  S  is  unitary,  then  0  is  a  J-unitary  matrix  [26],  [25],  [21]: 


0H  JO  =  j 


1  0 
0  -1 


Thus,  the  lossless  2-ports  provide  excellent  examples  of  unitary  and  J-unitary  ma¬ 
trices.  The  mappings  between  the  scattering  matrices  and  the  chain  matrices  are 
provided  by  Hasler  and  Neirynck  [18]: 


S  s21x 


—  detfS1]  sn 
— s22  1 


0  h_ ►  ^22 


6 12  det[@] 

1  —  9  21 


S. 


(4) 


10 


h 


12 


Figure  2:  Chain  of  2-ports  has  chain  scattering  matrix  0  =  @i@2- 


Although  the  2-port  has  a  scattering  matrix,  it  admits  a  chain  scattering  matrix  only 
if  .S'2 1  is  invertible.  These  notions  generalize  to  iV-ports.  A  matrix  0  is  J-unitary 
when 


QHje  =  J:  = 


Im  0 
0  — In 


The  collection  of  all  such  J-unitary  matrices  is  denoted  as  U(M,N).  Equation  4 
generalizes  to  the  map  S  :  U(M,N )  — >  U(M  +  N )  that  takes  a  chain  scattering 
matrix  to  its  corresponding  scattering  matrix  [21],  [8]: 


s|e] 


©120221  @11  —  0120221@21 
©221  —  0  221@21 


(5) 


The  matrix  map  S  is  well-defined  because 


[  Im 

1 

O 

'  0fi©n 

-  ©21021 

©fl©12 

—  ©21  ©22 

1 

O 

—In 

.  0f20n 

@22©21 

©f2©12 

@22©22 

0HJ0 


forces  both  @n  >  I m  and  ©22  >  In-  The  matrix  map  S  turns  the  CSD  and  unitary 
dilations  into  a  hyperbolic  CSD  and  J-unitary  dilations. 


7  Cosh-Sinh  Decomposition 


The  canonical  example  of  a  hyperbolic  matrix  is  given  by  Mendes  and  Ruas  [26] : 


H 


cosh(4/?)  sinh(,0) 
sinh('0)  cosh^) 


6W(1,1). 


This  example  is  a  special  case  of  the  HCSD. 
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Corollary  4  (HCSD)  Let  H  G  U(M,N)  with  M  <  N.  Then  there  are  unitary 
matrices  Tn,  Ti  G  U(M)  and  U22,  V22  G  U(N)  such  that 


H 


Tn  0 
0  U22 


cosh(T) 

sinh(T)  0 

sinh(T) 

cosh(T)  0 

0 

1 

5= 

1 

fe; 

O 

'  Ti 

0  ' 

0 

T22  _ 

H 


where  T  =  diag(Vh,  •  •  • ,  iPm)  for  ifm  >  0. 


Proof:  If  H  is  J-unitary,  Equation  5  makes 


W  =  S[H } 


H\  2  H22  h11-h12ht21h21 
I  Iff  ~Hf21H21 


unitary.  To  preserve  the  block  structure,  apply  the  CSD  as  follows: 


IT,  2  IT,, 
IT22  IT2 , 


Tn  0 
0  U22 


cos(0) 

—  sin(@) 

0 

sin(0) 

cos(0) 

0 

0 

0 

In-m 

'  Tn 

0  ' 

0 

T22  _ 

H 


where  0  =  diag(0, ,  •  •  • ,  9m)  for  0  <  6rn  <  ir/2.  Substitution  back  into  H  gives 


H 


Tn  0 

0  —  T22 


cos(0)  1 

tan(0)  0 

tan(@) 

cos(0)_1  0 

0 

1 

O 

Tn 

0 


0 


H 


-T22 


Relabeling  produces  the  hyperbolic  CS  decomposition.  /// 


Just  as  the  CSD  parameterizes  the  unitary  dilations,  the  HCSD  parameterizes  the 
hyperbolic  dilations. 


8  Parameterizing  Hyperbolic  Dilations 

A  hyperbolic  dilation  can  be  defined  in  several  ways.  One  approach  starts  with  matrix 
A  G  CMxM  such  that  AH A  >  Im ■  The  problem  is  to  find  matrices  B,  C,  and  D  such 
that 

Ha  [CD 

is  J-unitary.  If  N  >  M  and  A  has  the  SVD 

A  =  Tn  cosh(T)Tf, 
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then  the  HCSD  and  Corollary  1  arguments  give  that  all  hyperbolic  dilations  of  A 
with  D  ~  N  x  N  have  the  form 


A  B 
C  D 


Uu  0 
0  U22 


cosh  (Hi) 

skill  (Hi)  0 

skill  (Hi) 

cosh  (Hi)  0 

0 

l 

0 

'  kj , 

0  ' 

0 

U22  _ 

H 


where  U22,  V22  G  U{N). 
If  A  has  the  SVD 


A 


Uu 


II  0 

0  cosh  (Hi) 


H 


Vu 


with  0  <  <  •  •  •  <  i/)k,  then  Corollary  2  arguments  ensure  that  the  smallest 

hyperbolic  dilations  of  A  have  the  form 


'  A 

B  ' 

'  Uu 

0  ‘ 

C 

D 

0 

U22  _ 

IL 

0 

0 

0 

cosh  (H/) 

skill  (T) 

0 

skill  (Hi) 

cosh  (Hi) 

V11 

0 


0 


H 


V22 


where  U2 2,  V22  G  U{K). 

These  hyperbolic  dilations  are  obtained  by  mapping  the  hyperbolic  matrix  back  to 
a  unitary  matrix,  dilating,  and  mapping  back  to  a  hyperbolic  dilation.  Moreover,  just 
as  the  CSD  generalizes  to  the  hyperbolic  CSD,  other  decompositions  of  linear  algebra 
admit  hyperbolic  counterparts.  Array  processing  problems  led  to  a  hyperbolic  SVD 
[29],  a  hyperbolic  URV  decomposition  [35],  and  a  hyperbolic  approach  to  Kalman 
filtering  [16] .  The  general  principle  is  that  the  dilations  and  decompositions  obtained 
for  one  class  of  matrices  can  map  into  another  class  of  matrices. 


Implementing  these  decompositions  in  VLSI  has  been  a  research  topic  at  the 
Institute  for  Network  Theory  and  Circuit  Design,  Technical  University,  Munchen, 
Germany  [11].  Diepold  and  Pauli  [9]  started  with  the  Schur  Decomposition  of  in¬ 
definite  matrices.  They  realized  this  decomposition  was  part  of  the  more  general 
problem  of  embedding  a  passive  matrix  in  a  lossless  matrix  [10].  They  subsequently 
found  a  group-theoretic  approach  that  organizes  this  embedding  problem  in  a  signal¬ 
processing  context  that  admits  hardware  solutions  [12]. 


9  Antenna-Matching  Applications 

Figure  3  shows  a  180°  hybrid  chained  to  a  double  ferrite  antenna  and  loaded  with  a 
lossless  2-port.  The  hybrid  and  the  antenna  are  4-ports.  The  lossless  2-port  is  the 
designable  part  of  this  system.  The  ports  are  connected  by  lines  that  represent  the 
two  wires  that  attach  to  the  two  terminals  that  constitute  a  port. 
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Figure  3:  Hybrid-Antenna  chain  loaded  with  lossless  2-port  (solid). 


When  the  lossless  2-port  terminates  the  antenna,  the  resulting  hybrid-antenna 
4-port  is  converted  to  a  2-port  with  scattering  matrix 


ST  = 


•Si  1  S12 

S21  S22 


The  design  goal  is  find  a  lossless  2-port  that  forces  Su  and  s  \ 2  to  be  small.  The  pa- 
rameterizations  of  the  lossless  2-port  provide  performance  bounds.  For  this  particular 
antenna  under  consideration,  these  performance  bounds  show  that  this  2-port  loading 
cannot  simultaneously  force  Sn  and  s  1 2  to  be  small.  Rather  than  waste  time  trying 
to  load  this  antenna,  the  engineer  should  look  for  other  antennas  that  are  amenable 
to  2-port  loading. 


The  hybrid  has  ideal  scattering  matrix  [31,  Equation  7.101]: 


SH,  0  =  -- U 


3_ 

v/2 


0 

1 

1 

0 


11  O' 
0  0-1 
0  0  1 
-1  1  0 


However,  Sh,q  implies  that  the  ports  are  indexed  as  shown  in  the  right  side  of  Figure  4. 
To  chain  the  hybrid  to  the  antenna,  the  ports  are  numbered  using  the  hybrid  on  the 
left  side  of  Figure  4. 


The  renumbered  hybrid  has  scattering  matrix  [24,  Equation  1]: 


Sh  = 


J 

V2 


0 

0 

1 

1 


0 

0 

-1 

1 


1  1 ' 
-1  1 
0  0 
0  0 


This  scattering  matrix  Sh  has  chain  matrix: 

'll  0  O' 

=  _J_  -110  0 

H  y/2  0  0  -1  -1 

0  0  1-1 
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Figure  4:  Hybrid  port  numbering. 


For  reference,  the  mapping  between  the  scattering  matrix  S  and  its  chain  scattering 
matrix  0  is  provided  by  Kimura  [25,  Equations  4.11,  4.12]: 

q  _  ©120221  ©11  —  ©120221021 

©221  — ©22  — 1021 

and  [25,  Equation  4.5]: 

c  c  o— i  c  c  o— i 
^12  —  ^11^21  *^22  ^11^21 
Q—l  Q  Q— 1 

—  021  *^22  *^21 

Figure  5  plots  the  scattering  matrix  Sa  of  the  Double  Ferrite  Antenna.  Each 
element  of  the  antenna’s  scattering  matrix 

SA,ll  SA,2\  SA,31  SA,A\ 

C  SA,  12  Sa,22  $A,32  S  A  A2 

&A  = 

SA,13  Sa,  23  Sa,33  S.4,43 

.  Sa,14  Sa,24  Sa,34  S.4,44 

is  a  complex- valued  function  of  frequency  s A.mnij^ f  )  for  2  <  /  <  10  MHz  in  the 
complex  unit  disk. 

Although  the  scattering  matrix  Sa  is  a  4  x  4  matrix,  Sa  is  both  symmetric, 

SA  =  Si. 

and  centro-symmetric  (symmetric  across  the  cross-diagonal), 

0  0  0  1' 

0  0  10 
0  10  0' 

1  0  0  0. 

In  addition,  Sa  is  constant  along  the  diagonal  and  the  cross-diagonal.  Consequently, 
Sa  can  have  no  more  than  four  independent  functions. 
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Double  Ferrite  Antenna:  4x  4  Scattering  Matrix 


Figure  5:  Scattering  matrix  of  Double  Ferrite  Antenna. 


If  0^4  denotes  the  chain  scattering  matrix  of  the  antenna,  the  hybrid-antenna 
chain  has  scattering  matrix  computed  from  the  product  of  the  chain  matrices: 


S 


0n©n- 


Figure  6  plots  the  16  complex- valued  functions.  The  hybrid  is  destroying  some  of  the 
antenna’s  symmetry. 

Partition  the  scattering  matrix  of  the  4-port  into  2x2  blocks: 

a  Sn  S12 


Let  the  4-port  be  terminated  in  the  lossless  2-port  with  2x2  scattering  matrix  Sl- 
The  resulting  2-port  St  is  obtained  by  looking  into  Ports  1  and  2  of  the  Hybrid  while 
Ports  3  and  4  of  the  antenna  are  terminated  in  the  lossless  2-port  Sl' 

St  =  Su  +  S12SL(I  -  S22SL)~1S21  = 


Sll  Sl2 
S‘21  S  22 
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Hybrid/ Antenna  Cham:  4x  4  Scattering  Matrix 

1 - 1 - i - i —  - 


0.8 

0.6 

0.4 


-0.8 


-1 

-1  -0.5  0  0.5  1 

m  .  2<f<  1 0  MHz 

Figure  6:  Scattering  matrix  of  Hybrid-Antenna  chain. 


Observe  that  S  is  really  a  function  of  the  lossless  load  and  frequency: 


S(SL-ju) 


su(SL-,j2-Kf )  s12(SL]j2Tcf) 
s2i(SL;j2TTf)  s22(SL\j2nf) 


The  design  goal  is  to  simultaneously  minimize  sn  and  S21  •  One  multi-objective  func¬ 
tion  is  the  worst  performance  at  any  frequency: 

max{|sii(Sz,;  j2nf)\  :  2  <  f  <  10} 
max{\s2l(SL-,j2/Kf)\  :  2  <  /  <  10} 

If  U  denotes  a  class  of  lossless  available  to  the  designer,  the  multiobjective  optimiza¬ 
tion  problem  is 

min{7(5'L)  :  SL  £  IA}. 

A  lower  bound  on  this  performance  can  be  obtained  by  fixing  a  frequency: 


7  (SL)  = 


MSL)l 


l(SL-fo) 


Sn(SL]j2TTfo)\ 

s2i(SL]j2‘Kfo)\ 


<  7 (SL). 
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Although  S'l(j27t/o)  is  a  unitary  matrix,  not  all  unitary  matrices  may  be  parameter¬ 
ized  by  this  Sl-  That  is, 


{SL (j27r/o)  :  SLeU}CU( 2). 

This  set  inclusion  forces  the  inequality: 

min{7(S,  f0):  S  eU(  2)}  <  ”/(SL]  /0)  <  'j(Sl)- 

Consequently,  a  lower  bound  on  the  matching  performance  may  be  obtained  by  fixing 
the  frequency  and  sweeping  over  the  2x2  unitary  matrices  of  U{ 2).  These  unitary 
matrices  can  be  parameterized  by  Corollary  1 — any  S  e  W(2)  has  the  representation 


’  ej4> 11  0  ’ 

cos (9) 

—  sin(6l) 

‘  eA-  n  0  " 

0  e-^22 

sin(0) 

cos(0) 

0  e.A’22 

Figure  7  plots  7 (S',  /o)  over  a  dense  sampling  of  S'  G  7/(2).  The  plot  shows  that  it  is 
impossible  to  make  sn  and  ,s 2 1  simultaneously  small  with  lossless  loading  at  single  fre¬ 
quency.  The  performance  can  only  be  worse  over  the  frequency  band.  Consequently, 
Figure  7  tells  the  antenna  engineer  not  to  waste  time  with  any  lossless  loading  design 
for  this  antenna.  Rather,  the  antenna  must  be  redesigned  for  better  performance. 


Hybrid+ Antenna  sweeping  over  all  2x2  unitary  loads 


Figure  7:  Lower  performance  bound  of  lossless  loading. 
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